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Abstract. A theoretical investigation of the pulsation behavior of so-named "anomalous" Cepheids is presented. 
The study is based on nonlinear convective pulsation models with Z=0.0001 and 0.0004, mass in the range 1.3- 
2.2Mq and various luminosity levels. Based on these computations, we derive period, bolometric light curves 
and the edges of the instability strip, showing that a variation of the metal abundance from Z=0.0001 to 0.0004 
has quite small effects on these quantities. Then, using bolometric corrections and color-temperature transfor- 
mations, we are able to provide the predicted relations connecting pulsational properties (periods, amplitudes) 
with magnitudes and colors in the various photometric bands. The theoretical pulsational scenario is compared 
to observed anomalous Cepheids in dwarf spheroidal galaxies and, in particular, the predicted mass-dependent 
Period-Magnitude- Amplitude and Period-Magnitude-Color relations are used to estimate individual mass values, 
as well as to discriminate between fundamental (F) and first-overtone (FO) pulsators. 



1. Introduction 
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Oh, Among the various types of radial pulsating stars, the so-named "anomalous" Cepheids (AGs) are, perhaps, the least 
Q ' explored ones. They are rare in Galactic Globular Glusters, whereas they have been observed in all the nearby dwarf 
|H ' spheroidal galaxies (see e.g., Nemec et al. 1988, Nemec, Nemec & Lutz 1994, Mateo 1998, Pritzl et al. 2002, Dall'Ora 
^ et al. 2003, Baldacci et al. 2003) in which they have been searched for. These variables are intrinsically brighter than 
. . RR Lyrae stars and, when their apparent magnitude is scaled to the magnitude level of RR Lyrae stars observed in 
^ ■ the same system, they arrange along a strip with the less luminous ones (~ 0.5 mag brighter than RR Lyrae stars) 
^ ] at periods ^ 0.30 days and the more luminous ones (~ 2.0 mag brighter than RR Lyrae stars) at periods ~ 2.0 
5_] days. If compared with other luminous metal-poor variables such as Population II Gepheids (see the recent review by 
5^ : Wallerstein 2002), they appear significantly brighter, at fixed period. 

A general consensus exists in the literature (see, e.g., Bono et al. 1997a [B97], Gaputo 1998 and references therein) on 
the fact that AG variables, like RR Lyrae stars and Population II Gepheids, are metal-poor He-burning stars in the post 
Zero Age Horizontal Branch (ZAHB) evolutionary phase. As reviewed by Gaputo (1998), for low metal abundances 
{Z < 0.0004) and relatively young ages (< 4 Gyr) the effective temperature of ZAHB models reaches a minimum 
(logTe ^ 3.76) for a mass of about 1.0-1.2 Mq, while if the mass increases above this value, both the luminosity and 
the effective temperature start increasing, forming the so called "ZAHB turnover" . For larger metallicities, the more 
massive ZAHB structures have brighter luminosities but effective temperatures rather close to the minimum effective 
temperature. Within such an evolutionary scenario, AG variables appear to belong to the post-turnover portion of the 
ZAHB which crosses the instability strip at luminosity higher than RR Lyrae pulsators. Their origin is still debated and 
the most widely accepted interpretations are: (1) they are young (< 5 Gyr) single stars due to recent star formation; 
(2) they formed from mass transfer in binary systems as old as the other stars in the same stellar system. 

The general properties of observed AGs have been analyzed by Nemec, Nemec & Lutz (1994) and by Pritzl et al. 
(2002) who provided empirical Period-Luminosity (PL) relations in the B and V bands for both fundamental (F) and 
first overtone (FO) pulsators. However, as shown by the quoted authors, the mode assignment for these variables is 
not an easy task and the behavior in the PL, Period- Amplitude (PA) and Period-Golor (PG) planes does not always 
allow to discriminate between F and FO pulsation modes. 
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Table 1. Input parameters for the grid of AC pulsation models with y=0.24 and Z=0.0001 and 0.0004. 



M/Mq 



log L/Lq 



1.3 
1.6 
1.9 
2.0 

2.2 



1.82, 1.92 
1.88, 1.98, 2.08 
1.98, 2.08, 2.18 
2.08, 2.18, 2.28 

2.08, 2.28 



Based on previous nonlinear convective pulsating models and evolutionary tracks, B97 constrained the blue absolute 

magnitude and the period of metal-poor central He-burning pulsators more massive than RR Lyrae stars, finding a 
minimum mass for the occurrence of this kind of variables around ~ l.SM© with Z=0.0001 and ~ I.SMq with 
Z=0.0004. The theoretical limits of the distribution of these high mass pulsators in the AMs vs P plane, where 
AMb is the calculated Mb scaled to the ZAHB magnitude at the mean effective temperature of the RR Lyrae strip, 
conformed well to the observed data of most AGs and suggested a minimum mass of about 1.5 Mq. According to the 
position in the PL and PA planes, B97 confirmed that observed AGs are a mix of F and FO pulsators and derived two 
distinct non parallel B band PL relations for the two pulsation modes. These relations were significantly different from 
the empirical parallel ones derived by NNL, whereas they better conform with the analysis by Nemec et al. (1988) 
and Pritzl et al. (2002). 

However, B97 pulsation models did not include the most recent input physics and covered a restricted range of 
stellar parameters. For this reason, we computed new models including updated opacity evaluations and spanning a 
wider range of stellar parameters. In this paper we present the results of these computations and some preliminary 
tests of the predictive capability of the updated pulsational scenario. 

The organization of the paper is the following: in Sect. 2 we present the models, the results concerning the topology 
of the instability strip and the light curves. The transformation from theoretical to observational quantities is presented 
in Sect. 3, together with a short discussion of the different types of mean magnitudes and colors. In the same section, 
we give the relations connecting pulsational parameters (periods, amplitudes) to absolute magnitudes and colors. The 
comparison with observed AGs in dwarf spheroidal galaxies is briefly discussed in Sect. 4 and some concluding remarks 
close the paper. 

2. The pulsation models 

2.1. Periods and instability strip 

The basic physics underlying radial pulsation, that is the conversion of thermal energy into kinetic energy, is quite well 
known (e.g., see the nice review by Madore & Frcedman 1991) and, as early shown by Sandage (1958), the natural 
consequence of the Ritter's period-mean density relation and of the Stefan's law connecting luminosity, radius and 
effective temperature is that the pulsation period is a function of the luminosity, mass and effective temperature of the 
pulsator. Present models are computed using a nonlinear, nonlocal and time-dependent convective hydrodynamical 
code which has been already discussed in a series of previous papers (see, e.g.. Bono & Stellingwerf 1994; Bono et al. 
1997a, b) and is not repeated here. It seems of importance only to note that, with respect to linear envelope models, 
the nonlinear approach provides many more predictions beside the pulsation period, such as the luminosity variation 
along the pulsation cycle, namely the bolometric light curve;, and its morphology. Moreover, the inclusion of convection, 
which tends to quench pulsation, is fundamental to gain reliable information on both the blue and red edges of the 
pulsation region. Even though we are aware that further computations are needed to refine our knowledge of the 
pulsation phenomenon, our hydrodynamical code and physical assumptions have already shown a nice agreement with 
several observed properties of classical Gepheid (see Bono, Marconi & Stellingwerf 1999; Bono, Gastellani & Marconi 
2000b, 2002; Bono et al. 1999, 2002; Gaputo et al. 2000, 2001; Fiorentino et al 2002) and RR Lyrae models (see 
Bono, Gastellani. Marconi 2000a, Bono et al. 1997b, 2003, Marconi et al. 2003), and we are confident that they can be 
extended to other types of pulsating stars in order to construct a homogeneous pulsational scenario. It is also important 
to consider that the study of complex stellar systems, where a variety of pulsating stars can be observed, requires that 
the same input physics should be adopted for both pulsation and evolutionary computations. For these reasons, as for 
RR Lyrae and classical Gepheid models, we adopt the updated compilations of the opacity tables (Iglesias & Rogers 
1996) used in modern star model computations. Finally, we have considered two chemical compositions representative 
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of the metal-poor abundances of most of the observed AGs (Z=0.0001 and 0.0004, with ¥=0.24), a mixing-length 
parameter l/Hp=1.5, and the mass and luminosity levels suggested by He-burning evolutionary models (see Castellani 
& Degl'Innoccnti 1995), as reported in Tabic 1. 

For each mass and luminosity, the computations are performed by decreasing the effective temperature Tg by steps 
of 100 K, and the period from the bluest to the reddest pulsating model is derived^. Note that increasing (decreasing) 
by 100 K the effective temperature of the bluest (reddest) pulsators yields non-pulsating structures. In analogy with 
other well studied pulsators such as RR Lyrae stars (see van Albada & Baker 1971, Bono et al. 1997b, Marconi et 
al. 2003) and classical Cepheids (see Bono, Castellani & Marconi 2000b), linear regression through the results allows 
us to derive analytical relations connecting the period of both fundamental (F) and first-overtone (FO) modes to the 
intrinsic stellar parameters, namely mass, luminosity, and effective temperature. We obtain 

logPp = 10.879(±0.001) + 0.8181ogi-0.6161ogM- 3.309 logTe (la) 



logPFO = 10.198(±0.001) + 0.7751ogL- 0.5301ogM-3.1431ogre (16) 

where mass M and luminosity L are in solar units. As a whole, the period of FO pulsators can be fundamentalised as 
logPp^log Pfo+0.13, in agreement with the procedure currently adopted for RR Lyrae stars. 

For fixed metal abundance, mass and luminosity, FO pulsators are generally bluer than F pulsators, so that the 
blue edge of the first-overtone pulsation region (FOBE) and the red edge of the fundamental one (FRE) can be taken 
as representative of the boundaries of the whole instability strip, at least in the studied range of mass and luminosity. 
Giving the adopted steps of 100 K in our computations, we adopt as FOBE and FRE the effective temperature of 
the bluest FO model and of the reddest F model, increased and decreased by 50 K, respectively, with an intrinsic 
uncertainty of ± 50 K. As for the internal topology of the pulsation region, the computations show that the blue 
edge for fundamental pulsation (FBE) is bluer than the first-overtone red edge (FORE), suggesting the occurrence 
of a middle region of the instability strip where both the pulsation modes can be stable, similarly to what found for 
RR Lyrae stars (see van Albada & Baker 1973, Marconi et al. 2003 and references therein). Moreover, we find that 
the effective temperature of reddest FO models, at fixed mass and luminosity, appears to be sensitive to the location 
of the pulsation node within the convective stellar envelope, deserving a detailed analysis that will be presented in a 
forthcoming paper. 

The lines drawn in Fig. 1 show, left to right, the location in the HR diagram of the bluest FO, the bluest F and the 

reddest F pulsators, respectively, for the labelled mass and varying the luminosity imdcr the assumption 1/ Hp=l.5. 
We find that, at constant mass and luminosity, a variation of Z from 0.0001 to 0.0004 has a marginal influence on the 
instability strip topology, the only evident effect being the narrowing of the region where only FO pulsation is efficient, 
that is the region between the bluest FO and F models, when passing from Z=0.0001 to Z=0.0004. As a whole, in 
spite of the different opacity tables, the predicted topology of the instability strip is similar to the one found by B97 
for a stellar mass 1.5 M© and Z=0.0001 (see Fig. 1 in B97). 

Inspection of the results plotted in Fig. 1 shows that the entire pulsation region tends to move towards the blue if 
increasing the mass (at constant luminosity) or decreasing the luminosity (at constant mass). A linear regression to 
the results yields that, in the range Z=0.0001 to 0.0004, the limits of the AG instability strip can be approximated as 

logTe{FOBE) = 3.947(±0.006) - 0.052 log i + 0.042 log M (2) 



log Te (Fi?£) =3.898(±0.006)- 0.065 log L-F 0.058 log M (3) 

where the intrinsic uncertainty of ± 50 K on the FOBE and FRE temperatures is included. Concerning the fundamental 
blue edge (FBE), our models suggest logTe(FBE)-logTe(FRE) ~ 0.055± 0.005, for fixed mass and luminosity. Note 
that, as already discussed for RR Lyrae stars (see Caputo et al. 2000, Marconi et al. 2003), predictions concerning 
the FOBE, that is the boundary between non pulsating and pulsating stars at the hotter effective temperatures, are 
particularly solid. In fact, the uncertainties connected to the treatment of cx)nvcction are expected to be smaller in the 
blue region (high effective temperatures) of the instability strip and specific tests show that a variation of the mixing 
length parameter from l/Hp from 1.5 (the adopted value in this investigation) to 2.0 implies that the FOBE becomes 
redder by not more than 100 K whereas the FRE becomes bluer by at least 200 K. 



^ The entire set of results is available upon request to the authors. 
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2.2. Light curves and pulsation amplitudes 

For each investigated model, the variation of relevant parameters (luminosity, radius, radial velocity, effective temper- 
ature and gravity) along a pulsation cycle was provided. A sub-sample of bolometric light curves for Z=0.0001 (with 
M/M0=2.O and logL/L0=2.O8) is shown in Figs. 2a, b for fundamental and first overtone pulsators. respectively. The 
behavior of the curves shows that, at variance with the case of RR Lyrae stars, no clear information on the pulsation 
mode can be drawn from the light curve morphology. Indeed, even if the raising branch is steeper in the case of 
the fundamental mode, the light curves for both modes present a variety of morphological features, deviating from a 
sinusoidal shape but for the few hottest first overtone pulsators. 

Figure 3 shows bolometric amplitudes versus periods for selected models with Z=0.0001. In particular, the effects 
of varying the luminosity level at fixed stellar mass {2.2Mq. top panel) and of varying the mass at fixed luminosity 
(logL/L0=2.O8, bottom panel) are shown. As a whole, fundamental pulsators follow the well known behavior of 
fundamental RR Lyrae stars, with the amplitude increasing as the period decreases, for fixed mass and luminosity. 
Also the dependence on mass and luminosity is the same as found for RR Lyrae stars (sec Bono ct al. 1997b, Marconi 
et al. 2003): at fixed period, the amplitude increases with increasing the luminosity or decreasing the mass. Moreover, 
the rather linear dependence of the amplitude with period will allow us (see Section 3) to derive fundamental Period- 
Amplitude relations in analogy with RR Lyrae stars. 

As for first overtone pulsators, let us recall that, for any fixed mass and luminosity, a model with an effective 
temperature 100 K higher than the hottest pulsating model is not pulsating. This means that within S logP ~ 0.02- 

0. 03 to the left of the hottest (shortest period) FO model plotted in Fig. 3 the amplitude should reach vanishing values. 
On this basis, the amplitudes plotted in Fig. 3 appear to show the characteristic bell shape typical of first overtone 
RR Lyrae stars. 

However, in spite of these similarities, two main differences appear to exist between the behaviors of AGs and RR 

Lyrae stars in the Bailey diagram: 

1. in the case of AGs the maximum amplitudes reached by first overtone pulsators may be comparable with those of 
fundamental ones, whereas fundamental RR Lyrae pulsators always attain significantly higher amplitudes; 

2. the maximum of the bell defined by first overtone pulsators decreases when the stellar mass increases or the 
luminosity level decreases, following an opposite trend with respect to first overtone RR Lyrae stars (see Bono et 
al. 1997b for details). 

Gonsidering that observed AGs, even in the same stellar system, may have different mass and luminosity, with 
the more massive ones at brighter luminosities, one has that the discrimination of the AG pulsation mode from the 
analysis of the period-amplitude diagram might be not an easy task. As a fact, inspection of Fig. 3 reveals that the 
period-amplitude distribution of the more massive (2.2M0) and brighter (logL/L0=2.28) FO pulsators (dotted line 
connecting open circles in the top panel), is quite close to that of less massive (I.6M0) and fainter (logL/i0=2.O8) F 
pulsators (solid line connecting dots in the bottom panel). 

3. From theoretical to observational quantities 

In order to compare theoretical results with observations, the bolometric light curves have to be transformed into the 
photometric bands UBVRIJK. As in previous studies performed by our group, we use bolometric corrections and 
temperature-color relations provided by Gastelli, Gratton & Kurucz (1997a, b). 

Before proceeding, it is worth noticing that for observed pulsating stars we can only measure magnitudes and 
colors that are time-averaged quantities over a pulsation period, rather than the static values the stars would have 
were they not pulsating. For this reason, in consideration that observed variables are fitted with evolutionary tracks, 
it is important to estimate the mean values from pulsation models and compare them with the corresponding static 
quantities. A similar study has already been performed for RR Lyrae stars (Bono, Gaputo & Stellingwerf 1995, Marconi 
et al. 2003) and classical Gepheids (Gaputo, Marconi & Ripepi 1999) and we are now in the position to extend the 
analysis to AGs. 

With < Mi > and (Mj) indicating intensity- weighted and magnitude- weighted mean magnitudes, respectively, we 
show in Fig. 4 the difference between static My and mean magnitudes as a function of the pulsation visual amplitude 
Ay One has that intensity- weighted means are better representative of the static values than magnitude-averaged 
quantities, in both cases the discrepancy increasing toward higher amplitudes and less symmetric light curves. As 
far as the BV color is concerned. Fig. 5 shows that increasing the visual amplitude the magnitude-averaged mean 
colors {B — V) tends to be redder than the corresponding static values, whereas the opposite behavior holds for the 
intensity- weighted mean colors < B > — < V >. 

In general, the behaviors found for present models are in agreement with the results already found for RR Lyrae 
stars and classical Gepheids. In analogy with the procedure followed by Marconi et al. (2003) for RR Lyrae stars, we 
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derive the predicted relation to estimate static Vg magnitudes from observed magnitude and intensity averaged values, 
as given by 

Vs = 1.20 <V > -0.20(F) (4) 

with an uncertainty of ± 0.01 mag. This relation appears in substantial agreement with the results recently presented 
for RR Lyrae models with Z=0.001 (see Marconi et al. 2003). 

However, bearing in mind that the pulsation amplitude of models close to the boundaries (FOBE and FRE) of 
the instability strip approaches the zero value, the location of FOBE and FRE in the My-logP plane turns out to be 
quite independent of whether static or mean magnitudes are considered. As shown in Fig. 6, a linear regression to the 
models with l/Hp=1.5 yields the analytical relations 

Mv{FOBE) = -1.12 -2.68 log Pto - 1.80 logM (5) 

Mv{FRE) = -0.03 - 2.321ogPF - 1.92 logM (6) 

with cry=±0.06 mag, including the intrinsic ± 50 K uncertainty on FOBE and FRE effective temperatures. Adopting 

l/Hp=2.0, the zero-point in cq. (5) and (6) changes respectively by^ +0.1 mag and ^ —0.18 mag. As for the blue 
limit of fundamental pulsation (FBE), where the amplitude reaches quite large values, we consider intensity- weighted 
magnitudes, deriving that the predicted location of the FBE can be approximated as 

< Mv{FBE) >= -0.72 -2.61 log Pf - 1.71 logM (7) 

with ^^=±0.07 mag. 

3.1. Period-Magnitude-Color, Color-Color relations and Wesenheit functions 

The natural outcome of the period relation (cqs. la, lb) in the observational plane is the Period-Magnitude-Color 
relation (hereafter named PMC) in which the period is correlated with the pulsator absolute magnitude and color, for 
any given mass. The coefficients of the intensity-weighted PMC relations in the various bands are reported in Table 
2, while Fig. 7 shows the quite small scatter of the models around the predicted relations. 

According to these relations, one could estimate the mass range spanned by a sample of ACs at the same distance 
and with the same reddening with a formal uncertainty of ~ 2%, provided that periods and colors are well measured. 
If the distance and reddening are independently known, then the mass absolute values can be determined too. 

We have also investigated the relations among intensity-averagcid colors, from optical to near infrared bands. 
As given in Table 3, these relations have an intrinsic dispersion of the order of few thousands of magnitude, but 
unfortunately they appear to be of no use for what concerns the reddening determination. As a fact, apart from the 
< B > — < V > versus < V > — < R > case, the slope b of the predicted relations is very close to the reddening 
vector r, as estimated using a typical extinction law (see e.g. Cardelli, Clayton & Mathis 1989). 

In order to complete the pulsational framework, we have finally considered the reddening free Wasenheit functions 
(Madore 1982) and their dependence on periods. As already discussed in several papers on classical Cepheids (sec, 
e.g., Madore & Freedman 1991, Tanvir 1999, Caputo, Marconi & Musella 2000), the intrinsic variation due to the 
finite width of the instability strip is almost similar to the effects of interstellar extinction, and the color term in the 
Wasenheit functions [W{B,V) = My-3.1(P- V"), W{V,I) = My -2.54(^-7), etc.], thought to cancel the reddening 
effect, helps to reduce the dispersion of magnitudes at a given period. On this basis, the Period- Wesenheit relations 
(hereafter referred to as PW) is used to get reliable estimates of the intrinsic distance modulus. 

This effect is shown in Fig. 8, while the coefficients of the predicted intensity weighted PW relations, as derived 
by fundamental models and first-overtone models with fundamentalised periods, are reported in Table 4 for selected 
photometric bands. Note that the use of these reddening-free relations to get the mass range of a sample of ACs at the 
same distance yields a formal uncertainty of ~ 5%, provided that periods and colors are well measured. On the other 
hand, if the mass of the variable is independently known, these relations can provide its intrinsic distance modulus. 

3.2. Period-Magnitude-Amplitude relation 

Once the bolometric light curves are transformed into the observational plane, pulsation amplitudes in the various 
photometric bands arc derived. As already noticed in Sect. 2, the behavior of fundamental pulsators in the Period- 
Amplitude diagram is rather linear, but with a non negligible dependence on both the pulsator mass and luminosity. 
This occurrence suggests the possibility to derive a mass-dependent relation between period, absolute magnitude and 
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Table 2. Coefficients of the mass-dependent Period-Magnitude- Color relations. 



Mode a /3 7 (5 av 



<Mv>=a + /31ogP-F7[< B > - <V >]+ SlogM 
F -1.56 -2.85 3.51 -1.88 0.01 

FO -1.92 -2.90 3.43 -1.82 0.02 

< Mv >= a + /31ogP-|-7[< V > - < R>]+ SlogM 
F -1.69 -2.90 5.39 -1.89 0.01 
FO -2.09 -3.05 5.28 -1.75 0.01 

< Mv >=a + piogP+j[< V > - < I >] + SlogM 
F -1.83 -2.93 2.73 -1.89 0.01 
FO -2.24 -3.09 2.67 -1.73 0.01 

< Mv >= a + /31ogP+7[< V > - < J >] + SlogM 
F -1.95 -3.01 1.78 -1.87 0.01 
FO -2.38 -3.17 1.75 -1.70 0.01 

< Mv >=a + l3logP+^[< V > - < K >] + SlogM 
F -1.99 -3.05 1.31 -1.86 0.01 
FO -2.40 -3.22 1.29 -1.68 0.01 

< Mv >=a + /31ogP+7[< J > - < K >]+ SlogM 
F -1.97 -3.07 2.17 -1.85 0.01 
FO -2.41 -3.24 2.11 -1.66 0.01 



Table 3. Coefficients of intensity averaged Color-Color relations for fundamental models at Z=0.0001. The last column 
gives the reddening vector, as estimated using a typical extinction law (see text). 



Mode a bar 



<B> 


- <V >■- 


^a + b[<V > 


- <R>] 


F 


-0.027 


1.509 0.002 


2.03 


FO 


-0.016 


1.461 0.004 


2.03 


<B> 


-<V> 


= a + b[<V> 


-<7>] 


F 


-0.060 


0.755 0.002 


0.82 


FO 


-0.046 


0.726 0.004 


0.82 


<B> 


- < I" > 


= a + b[<V > 


- < J >] 


F 


-0.077 


0.478 0.004 


0.45 


FO 


-0.062 


0.461 0.007 


0.45 


<B> 


-<V>-- 


= a + b[<V> 


-<K> 


F 


-0.070 


0.347 0.005 


0.36 


FO 


-0.056 


0.335 0.008 


0.36 



amplitude (PMA). A linear regression through the fundamental models with Z=0.0001 and 0.0004 gives the relation 
(see also Fig. 9) 



logPp* = logPp -I- 0.41 <Mv> -1-0.77 log M = 
0.01(±0.05)-0.188Ay 



(8) 
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Q /3 7 av 



< W{B,V) >= 


Mv > 


-3.10[< B > - 


<V >] 


-1.43 -2.77 


-1.83 


0.04 




< W{V,I) >= 


:< Mv > 


-2.54[< V > - 


<I>] 


-1.74 -2.94 


-1.83 


0.04 




< W{V,K) >= 


:< Mv > 


-1.13[< V > ~ 


<K>] 


-1.73 -2.93 


-1.83 


0.04 





that again should allow to estimate the mass range spanned by fundamental AGs at the same distance and with well 
measured periods and amplitudes, but with a formal uncertainty of ^ 15% which is larger than that related with the 
PMC and PW relations. However, it is of interest to note that, since the PMA relation only holds for fundamental 
pulsators, a comparison with the results inferred from the PMC or PW relations should provide a way to discriminate 
the pulsation mode of observed ACs (see below) . 

4. Comparison with observations 

4.1. Anomalous Cepheids in dwarf spheroidal galaxies 

As clearly shown by the evolutionary models (see Castellani & Degl'Innocenti 1995, B97, and the recent Pisa 
(www.mporzio.astro.it/ marco/GIPSY/homegipsy.html) and Padova ( http: / /pleiadi.pd. astro. it) libraries), metal-poor 
central He-burning structures more massive than ^ I.SA^q evolve with a characteristic "blue loop" whose luminosity 
and effective temperature increases when increasing the mass. Since the effective temperature at the edges of the 
instability strip is a function of mass and luminosity, it follows that one can easily estimate whether a star with a 
given mass is expected to cross the pulsation region becoming a variable star. Moreover, since at fixed metal content 
the He-burning pulsators of a given mass would have a quite defined range of intrinsic luminosity, the mass-term in 
the above predicted relations could be removed to provide "evolutionary" relations. 

However, from the point of view of stellar evolution one should consider that the mass/luminosity relation of central 
He-burning structures is dependent on the metal content, as well as that, even in presence of stellar models computed 
with the most updated input physics, the He-burning luminosity of a given mass depends on the adopted efficiency of 
core overshooting. As shown in Fig. 10 for two selected mass values, canonical He-burning models computed with an 
inefficient overshooting (Pisa library) are brighter by ^ 0.25 mag with respect to overshooting models (Girardi et al. 
2000). However, if for the given mass the predicted FRE is taken into consideration, then the inclusion of overshooting 
seems to leave almost unvaried the minimum mass for the occurrence of AC variables. On this basis, in this paper we 
wish to present only the results of pulsation models, leaving mass and luminosity as free parameters and deserving to 
a forthcoming paper the connection to stellar evolution theory and the construction of synthetic relations based also 
on the evolutionary times. 

In this section, we are going to apply the theoretical relations presented in the previous section to observed AC 
samples in dwarf spheroidal galaxies with the purpose to estimate their mass. We searched in the literature for bona 
fide AC pulsators with well measured periods, magnitudes and amplitudes and we found the samples reported in 
Table 5 where the columns give the host galaxy, the number of measured ACs, the available observed quantities, the 
intrinsic distance modulus /iq, the E{B — V) reddening (both /io and E{B — V) as listed by Pritzl et al. 2002) and the 
corresponding references. 

A first comparison between theoretical results and observations is presented in Fig. 11 where observed ACs are 
plotted in the My-^ogP diagram together with the predicted boundaries (FOBE and FRE) of the instability strip at 
1.3Mq (dashed line) and 2.2Mq (solid line). We notice that the general behavior of observed pulsators is in substantial 
agreement with theory. Specifically, the comparison between the observed distribution and the predicted edges confirms 
that the brightest variables should have larger masses than the faintest ones. 

For almost all the dwarf galaxies reported in Table 5, except Carina and Fornax, a mode discrimination has been 
already attempted (see Pritzl et al. 2002), as listed in the first column of Table 6. Based on this classification and using 
the intrinsic distance moduli and reddening values listed in Table 5, we proceed in estimating the mass of individual 
variables from the predicted PMC relations given in Table 2. The results are given in columns (7) and (8) in Table 6. 
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Table 5. Selected parameters of anomalous Cepheids in dwarf spheroidal galaxies. 



host galaxy 


AC number 


observed quantities 


Mo 


E{B - V) 


source 


And VI 


6 


V,B,Av 


24.45 mag 


0.06 mag 


Pritzl et al. 2002 


Leo II 


4 


V,Av 


21.59 mag 


0.02 mag 


Pritzl et al. 2002 


Draco 


8 


V,B,Av 


19.49 mag 


0.03 mag 


Pritzl et al. 2002 


Carina 


8 


V,B,Av 


20.14 mag 


0.04 mag 


DaU'Ora ot al. 2003 


Sculptor 


3 


V,Av 


19.56 mag 


0.02 mag 


Kaluzny et al. 1995 


Sextans 


6 


V,B,Av 


19.74 mag 


0.03 mag 


Mateo et al. 1995 


Fornax 


17 


v,i 


20.70 mag 


0.03 mag 


Bersier & Wood 2002 



Furthermore, for F pulsators only, the mass is estimated also from the predicted PMA relation (sec cq. 8), yielding 
the values listed in column (9) in Tabic 6. It seems worth noticing that the stellar masses inferred in this paper appear 
in substantial agreement with the mass range of AGs, as predicted by evolutionary models of not-too-old metal-poor 
HB structures (see Castellani & Degl'Innocenti 1995, B97, and the quoted Pisa and Padova libraries). 

The mass values listed in Table 6 for the fundamental pulsators arc plotted in the upper panel in Fig. 12, together 
with the formal errors as given by the intrinsic dispersion of the predicted PMC and PMA relations. The same intrinsic 
dispersions yield that the two mass values should agree at a level of Slog M ~ ±0.07 (solid lines). Allowing for further 
uncertainties as due to photometric errors and light curve fitting, the results plotted in the upper panel in Fig. 12 
confirm that for all the fundamental candidates, except V157 in Draco (FO pulsator?), the mass values predicted by 
the two relations agree with each other within the errors, even though the mass suggested by the PMC relation appears 
systematically smaller (Alog M ^ —0.09) than the PMA results. This effect could be connected to the theoretical 
uncertainty on the I /Hp parameter which affects pulsation amplitudes and in turn the PMA relation. This also means 
that one could in principle calibrate this free parameter by requiring the agreement between the PMA results and 
those inferred by the PMC, which is independent of the adopted I /Hp value. Finally, it is of interest to notice that the 
difference between the two mass values is independent of errors on the apparent distance modulus, whereas it depends 
on the adopted reddening. Moreover, the mass absolute values listed in Table 6 should vary by about ± 12% with 
6iJ.=± 0.1 mag. 

Let us now assume that also the FO variables are actually fundamental pulsators. Using again the fundamental 
PMC and PMA relations, one derives the quite unrealistic large masses given in brackets in Table 6. Moreover, as shown 
in the lower panel in Fig. 12, the difference between the two mass values is now significantly larger (AlogA^f ~ —0.17) 
than the allowed uncertainty, with the exception of V9 in Sextans. In summary, observed data seems to support the 
hypothesis that, if the variable is a fundamental pulsator, then the mass values inferred from the predicted PMC and 
PMA relations arc internally consistent within the errors. However the location of V9 in Sextans, if this variable is a 
FO pulsator, would suggest that this condition is necessary but not sufficient to discriminate the fundamental mode. 
On this ground, according to the mass values listed in Table 6, let us adopt that V158 and V182 in Carina are bona 
fide fundamental pulsators, whereas V29, V187, V193 and V205 should pulsate in the first-overtone mode. As for V33 
and V14, no discrimination is presently made. 

Once the pulsation mode is discriminated, F and FO variables in Table 6 can be used to get some useful empirical 
relations. We show in Fig. 13 that the Period-Magnitude (PM) relations for fundamental (including V158 and V182 
in Carina) and first-overtone pulsators (including V29, V187 and V193 in Carina) can be linearly approximated as 

< Mv >F= -0.69(±0.14) - 2.731ogPF, (9a) 

< Mv >FO= -1.56(±0.25) - 2.951ogPFO, (96) 

which are in fair agreement with the relations recently provided by Pritzl et a. (2002). 

Moreover, using only the variables with BV data, we plot in Fig. 14 the absolute < WBV > quantities versus 
period. We derive that the PW relations are 

< WBV >F= -1.72(±0.20) - 3.34 log Pf, (10a) 
< WBV >FO= -2.26(±0.20) - 2.441ogPi;'o (106) 
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Finally, using only fundamental variables, we show in Fig. 15 that the empirical fundamental PMA relation can 
be written as 

<Mv>F= -0.37(±0.18) - 2.57(logPF + 0.188^y) (11) 

In these figures, the arrows depict the variables V157 in Draco (filled triangle) and V9 in Sextans (open dot), while 
the small crosses show V14 and V33 in Carina for which there is no clear classification. By inspection of the three 
figures, one derives that VI 57 in Draco follows the fundamental PM and PMA relations, but the PW relation for 
FO pulsators, likely suggesting that its measured BV color is too red. As for V9 in Sextans and the two variables in 
Carina, one derives that they follow more closely the FO relations rather than the fundamental ones. 

In conclusion, the mode discrimination for a sample of anomalous Cepheids requires the knowledge of periods, 
colors and amplitudes. Assuming that the variables listed in Table 6 are representative of the average properties of 
ACs, fundamental or first-overtone candidates should agree with all the above empirical relations. 

In order to show the importance of both color and amplitude measurements, let us take into consideration the ACs 
observed in Fornax by Bersier & Wood (2002). Using the VI data listed in Table 7, the predicted PMC relations for F 
and FO pulsators (see Table 2) yield the mass values listed in columns (5) and (6), respectively. As a whole, the results 
are consistent with the mass values listed in columns (6) and (7) in Table 6, except the variable J024002.7 which 
appears to be too massive. Excluding this variable, we show in the upper panel of Fig. 16 the comparison between 
observed magnitudes, corrected with a visual distance modulus /xy=20.79 mag (from the distance and reddening values 
reported in Table 5), and the above empirical PL relations for F and FO variables. One has that most of the variables 
appear to be fundamental pulsators. 

On the other hand, since our pulsation models give 

< WVI >= -G.22(±0.001) + 1.03 < WBV > (12) 
eqs. (10a) and (10b) can be easily transformed into 

< WVI >F= -1.93(±0.20) - 3.341ogPf , (13a) 



< WVI >FO= -2.52(±0.20) - 2.51 log Pto (135) 

giving us a further opportunity for determining the pulsation mode. As shown in the lower panel of Fig. 16, there are 
now at least six variables which fit quite well the empirical FO relation. Such a discrepancy may be due to poor light 
curves in the / band, but without the further constrains provided by the PMA relation, almost for half of the whole 
sample of variables in Fornax we are unable to clearly discriminate between F and FO pulsators. 

4.2. The case of Leo A 

Recent results of a search for short-period variable stars in the metal-poor {Z ~ 0.0004) Leo A dwarf irregular galaxy 
are presented by Dolphin et al. (2002). These authors have found 66 variables brighter than RR Lyrae stars, noticing 

that there is some ambiguity as to whether these objects should be classified. Eventually, they conclude that these 
variables could reproduce the natural extension to low metal abundance of classical Cepheids, rather than being 
"anomalous" Cepheids. 

We show in Fig. 17 the My-logP diagram of variables with unambiguous period and quality rating > 3, adopting 
an apparent distance modulus of 24.6 mag (see Dolphin et al. 2002). Filled and open symbols refer to F and FO 
pulsators, respectively, as given in Dolphin et al. (2002). The comparison with the data in Fig. 11 shows that the 
"faint" (My > —0.5 mag) variables observed in dwarf spheroidal galaxies are absent in Leo A. This evidence would 
suggest that the variables in Leo A are on average more massive than ACs in dwarf spheroidal galaxies. 

As a fact, the comparison with the predicted boundaries of the instability strip with M=2.2M0 (solid lines), 
discloses a quite reasonable agreement at the low luminosities, but the occurrence of bright (My ~ —1 to —1.5 mag) 
first overtone candidates located to the left of the predicted FOBE indicates the presence of even larger masses. In 
order to verify such a hypothesis, we have computed few models with 4.OM0, finding that the predicted FOBE at such 
a mass (dotted line) fits quite well the observed distribution. Based on such a result, we believe that the variables in Leo 
A are a mix of "anomalous" and "classical" Cepheids. A critical analysis of these variables requires both pulsational 
and evolutionary models, and this issue will be addressed in a forthcoming paper (Caputo et al. in preparation). 
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5. Concluding remarks 

Based on nonlinear, nonlocal and time-dependent convective pulsating models, we have presented a pulsational scenario 
for the analysis of "anomalous" Cepheid (AC) stars, showing that a variation of metallicity from Z=0.0001 to 0.0004 
has a very mild effect on the edges of the instability strip, as well as on the predicted relations connecting the pulsational 
(period, amplitude) to the intrinsic stellar parameters (mass, luminosity, effective temperature). The main results of 
the present study can be summarized a follows. 

1. We derive the pulsation relation connecting period with the star structural parameters M, L,Te. 

2. We give the predicted relations of the boundaries of the pulsation region in the logL-logTg and My-logP planes, 
for fixed mass. These relations, together with evolutionary models computed under different physical conditions 
(e.g., with or without overshooting) are the basic ingredients to evaluate whether a star is expected to evolve into 
the instability strip during its central He-burning phase. 

3. We show that the comparison between the mean magnitudes and colors, as based on theoretical light curves 
transformed into the observational bands, and the corresponding static values, follows the general trend also shown 
by other classes of pulsating stars (RR Lyrae, classical Cepheids). 

4. We derive the predicted Period-Magnitude-Color (PMC) relations, showing that it can be used to estimate, with a 
formal uncertainty of ~ 2%, the mass range spanned by a sample of well-measured variables at the same distance 
and with the same reddening. If the distance and reddening are independently known, then the mass absolute 
values can be determined too. 

5. As usually adopted for classical Cepheids, an important tool to derive absolute distances is given by the theoretical 
reddcning-frcc Wcsonhcit functions. These relations (PW) can also be used to get the mass range of a sample of 
ACs at the same distance with a formal accuracy of ~ 5 %. 

6. As already known for RR Lyrae stars, the amplitude of fundamental pulsators is linearly dependent on the pe- 
riod (logarithmic scale), for fixed mass and luminosity, providing a mass-dependent Period-Magnitude- Amplitude 
(PMA) relation which again can be adopted to estimate the mass range spanned by fundamental variables at the 
same distance, but with a larger formal uncertainty (~ 15%) with respect to the other methods. 

7. Since the PMA relation only holds for fundamental pulsators, whereas PMC and PW relations can be applied 
both to the fundamental and the first overtone mode, the comparison between the masses evaluated by adopting 
the different kinds of relations may help to attempt a discrimination of the pulsation mode. Moreover, we show 
that the difference between the mass values inferred from PMC and PMA relations is independent of errors on the 
distance modulus, whereas the mass absolute values should vary by about ± 12% with (5/x=± 0.1 mag. 

The comparison between the predicted edges of the instability strip and the observed distribution of ACs belonging 
to a number of dwarf spheroidal galaxies seems to support the predictive capabilities of pulsation models. The only 
exception is found for few first overtone variables in the irregular galaxy Leo A, that fall to the left of the predicted 
blue boundary reproducing the other variables, suggesting a mass of the order of 4 Mq. 

By considering dwarf spheroidal galaxies for which a mode discrimination of ACs is available in the literature and 
by applying the predicted PMC and PMA relations to these variables, we are able to conclude that if the variable is a 
fundamental pulsator, then both PMC and PMA yield mass values which are internally consistent within the errors. 
On this basis, we find that only two variables in Carina are bona fide fundamental pulsators, whereas the remaining 
variables should pulsate in the first-overtone mode. 

On the basis of the quoted mode discrimination empirical PL and PW relations for fundamental and first overtone 
ACs in dwarf galaxies, as well as the PA relation for fundamental variables, are derived. 

Finally, we notice that, as remarked by our referee, the intrinsic errors of the results listed above are 
valid in the context of the adopted treatment for the convective transfer. Indeed, even if our formalism 
has the advantages of being nonlocal and time dependent, nevertheless it has turned out recently 
that the inclusion of additional equations in the Reynolds momentum equations (see e.g. Canuto &: 
Dubovikov 1998) improves the treatment of turbulence. 
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Table 6. Estimated mass values for F and FO anomalous Cepheids, as inferred by PMC and PMA relations. The 
mode pulsation given in the first column is taken by previous studies, while that given in the last column is determined 
in the present paper (see text). 
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Table 7. Estimated mass values anomalous Cepheids in Fornax, as inferred by PMC relations. The mode pulsation 
given in the last column is determined in the present paper (see text). 
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Fig. 1. Dependence of the FOBE (dashed line), FBE (dotted line) and FRE (solid line) on mass and luminosity, for 
the two labelled metal abundances. 

Fig. 2. Theoretical bolomctric light curves for a subsample (M=2M0 logL/Z/0=2.O8) of fundamental (a) and first 
overtone (b) models with Z=0.0001. 

Fig. 3. Bolometric amplitudes versus periods for selected models with Z=0.0001 at varying the luminosity level at 
fixed stellar mass (2.2M0, top panel) and the mass at fixed luminosity (logi/Z/0=2.O8, bottom panel). 

Fig. 4. Difference between static My and either intensity weighted (top panels) or magnitude averaged (bottom panels) 
mean magnitudes as a function of the pulsation amplitude for both fundamental (F) and first overtone (FO) models. 

Fig. 5. The same but for the B-V color. 

Fig. 6. Limits of the pulsation region (FOBE and FRE) and blue limit for fundamental pulsation (FBE) in the 
My-logP plane. The periods are normalized to the same mass (I.9M0), according to the predicted dependence (see 
text). 

Fig. 7. PMC relations in the labelled photometric bands for both fundamental (filled circles) and first overtone(open 
circles) models. The solid line represents the linear regression through the data (see Table 2). 

Fig. 8. PW relations in the labelled photometric bands for fundamental (filled circles) and first overtone(open circles) 
models with fundamentalised periods. The solid line represents the linear regression through the data (see Table 4). 

Fig. 9. PMA relations for fundamental models. The solid line represents the linear regression through the data (see 
eq. 8). 

Fig. 10. Central He-burning evolutionary tracks for the labelled masses and metal content as computed without (solid 
line) or with overshooting (dashed line). The predicted FRE is also shown. 

Fig. 11. Comparison between the location of observed ACs in the My-logP diagram and the predicted boundaries of 
the instability strip at 1.3M0 (dashed line) and 2.2M0 (solid line). 

Fig. 12. (upper panel) - Comparison between the mass values of observed fundamental ACs (symbols as in Fig. 11), 
as inferred by the predicted PMA and PMC relations. The solid lines depict the allowed uncertainty, (lower panel) - 
The same, but for observed first-overtone ACs treated like fundamental pulsators. 

Fig. 13. Empirical My-logP relations for observed fundamental (upper panel) and first-overtone ACs (lower panel). 
The small arrows show V157 in Draco and V9 in Sextans discussed in the text. The solid lines are from eqs. (9a)-(9b) 
in this paper. 

Fig. 14. As in Fig. 13, but with empirical < WBV > quantities. The solid lines are from eqs. (lOa)-(lOb) in this 
paper. 

Fig. 15. Empirical PMA relation for fundamental AC. The small arrows show V157 in Draco and V9 in Sextans and 
the solid line is from eq. (11) in this paper. 

Fig. 16. (upper panel) - My versus logP for ACs in Fornax. The solid lines are the same drawn in Fig. 13. (lower 
panel) - < WVI > versus logP for the same variables. The solid lines are from eqs. (lOa)-(lOb) in this paper. 
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Fig. 17. My versus logP for F (filled dots) and FO (open dots) variables in Leo A with unambiguous period and 
quality rating > 3, adopting an apparent distance modulus of 24.6 mag (see Dolphin et al. 2002). The solid lines 
represent the predicted boundaries of the instability strip with M=2.2M0, while the dotted line is the predicted 
FOBE at M=4.OM0(see text for details). 
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